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Abstract 

Hamiltonian light-front field theory can be used to solve for hadron states 
' in QCD. To this end, a method has been developed for systematic renormal- 

ization of Hamiltonian light-front field theories, with the hope of applying 
the method to QCD. It assumed massless particles, so its immediate appli- 
cation to QCD is limited to gluon states or states where quark masses can 
\q ■ be neglected. This paper builds on the previous work by including particle 

masses non-perturbatively, which is necessary for a full treatment of QCD. 
We show that several subtle new issues are encountered when including 
■ masses non-perturbatively. The method with masses is algebraically and 

conceptually more difficult; however, we focus on how the methods differ. 
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We demonstrate the method using massive (j) 3 theory in 5+1 dimensions, 
which has important similarities to QCD. 



1 Introduction 



The use of a Hamiltonian light-front formalism may simplify the solution of quantum chromody- 
namics (QCD) by allowing us to make a convergent expansion of hadron states in free-particle 
Fock-space sectors. The Fock-space expansion will rapidly converge if the Hamiltonian satisfies 
certain conditions O. The Hamiltonian can then be used to solve for approximate hadron states. 

Inspired by the work of Dyson @ , Wilson || , Glazek and Wilson 0j , and Wegner H , significant 
work has been done to perturbatively derive light-front Hamiltonians that satisfy these conditions 
in the full Fock-space, neglecting zero modes [6-13]. The method developed by Allen and Perry 
[]l|] includes the scale dependence of the coupling and can be used to systematically renormalize 
light-front Hamiltonians, fixing all non-canonical operators, in principle to all orders. 

In this method the theory is regulated by placing on the Hamiltonian a smooth cutoff on change 
in free mass. The cutoff violates several physical principles, preventing renormalization exclusively 
through the redefinition of masses and canonical couplings. Renormalization must be completed by 
requiring the Hamiltonian to produce cutoff-independent physical quantities, and by requiring it to 
obey the physical principles of the theory that are not violated by the cutoff. These requirements 
completely fix the Hamiltonian so that it will give results consistent with all the physical principles 
of the theory, even those violated by the cutoff. The most powerful characteristic of this approach 
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is that it systematically "repairs" the theory and retains only the fundamental parameters of the 
theory. 

We generalize this method of renormalization to include particle masses and demonstrate it 
using massive 4> 3 theory in 5+1 dimensions. This theory is asymptotically free and its diagram 
structure is similar to QCD, which make it a good perturbative development ground. It is straight- 
forward to extend the method for massless theories developed in Ref. [|J to calculate QCD quantities 
for which particle masses are unimportant, such as the low- lying glueball spectrum [|14]|. In this 
paper, we show how to incorporate particle masses non-perturbatively as a necessary step toward 
a treatment of full QCD. 

2 Review of General Formalism 

In this section we introduce some of the notation developed in Ref. [jl[] and outline the method. 
Formalism that is necessary for a detailed understanding of this method but that we do not repeat 
in this paper can be found in this earlier work. This includes the use of a unitary transformation 
to determine how the Hamiltonian changes with the cutoff ||, the use of physical principles to 
restrict the form of the Hamiltonian, and the details of how to compute matrix elements of the 
Hamiltonian. 

We want to find the regulated invariant-mass operator, Ai 2 (g A ,m, A), which is trivially related 
to the Hamiltonian. It can be split into a free part (which contains implicit mass dependence) and 
an interacting part: 

M 2 (g A ,m,A) = Ml cc + M 2 nt {g A ,m,A). (1) 

Since the method treats Mf nt (g A , m, A) perturbatively, we put the particle-mass term in M 2 ree , to 
treat it non-perturbatively; however, Aif nt (g A ,m, A) will still have mass dependence. The matrix 
elements of M 2 (g A ,m, A) are written 

(F|^ 2 ( 5A ,m,A)|I) = {F\Ml cc \I) + (F\Ml t (g A ,m,A)\I) 

A 2 

= M 2 F (F\I) + e-^(F\V(g A ,m,A)\I), (2) 

where \F) and \I) are eigenstates of the free invariant-mass operator with eigenvalues Mp and M 2 , 
and Api is the difference of these eigenvalues. V(g A ,m, A) is the interacting part of the invariant- 
mass operator with the Gaussian cutoff factor removed and is called the "reduced interaction." 
With light masses Apj will be small in a limited part of phase space. This means there will be 
non-perturbative effects that must be dealt with through the diagonalization of the mass operator 
rather than its perturbative renormalization. 

We expand V(g A ,m, A) in powers of the running coupling, g A : 

oo 

V(g A ,m,A) = Y / 9 r A V^\m,A), (3) 

r=l 

where is the canonical interaction and the V^ r - 2 \m, A)'s are non-canonical interactions. These 
non-canonical operators can be thought of as counterterms in a traditional approach. Note that g A 
implicitly depends on m. We use a unitary transformation to relate A4 2 (g A ,m, A) to A4 2 (g A , , m, A'), 
where A' > A. This yields the relation 

r-1 

V (r) (m, A) - V (r) (m, A') = 5V ^(m, A, A') - ^ B r _ s<a V {r ' s \m, A) , (4) 

s=2 
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where 5V^ r \m, A, A') is the 0(g r r ) change in the reduced interaction and the -B r _ SjS 's are functions 
of m, A, and A' that contain information on the scale dependence of the coupling. Since the scale 
dependence of the reduced interaction comes from g A and the VM(m,A)'s [See Eq. (|)], Eq. (@) 
simply states that if we subtract from SV^ r \m, A, A') the contribution due to the scale dependence of 
the coupling, then we are left with the contribution due to the scale dependence of the V^ r \m, A)'s. 

If there is a part of V^ r \m, A) that is independent of the cutoff, it will cancel on the left- 
hand-side of Eq. (|j). For this reason, we split V^ r \m, A) into a part that depends on the cutoff, 
V^^(m, A), and a part that is independent of the cutoff, V^\m): 

V {r \m, A) = V^{m, A) + V$(m). (5) 

This division can be made with no ambiguity because we are assuming approximate transverse 
locality. Solving for both V^^m, A) and Vo^(m) is necessary to find the invariant-mass operator. 



3 Addition of Particle Masses 

In our renormalized scalar theory m is the physical particle mass to all orders in perturbation 
theory. In a confining theory m would be considered to be the particle mass in the zero-coupling 
limit. Since the mass is being treated non-perturbatively, it must be included in the free part of 
M 2 (g A ,m,A) in Eq. (|). This leads to an altered unitary transformation and fundamental changes 
in the renormalization procedure. 

The changes in the procedure are discussed in the next three subsections. The redefinition 
of the coupling (Sec. |3.l[) is straightforward. In Sections |3^ and |3.3| , we present the expressions 



for the matrix elements of V^^m, A) and V^^m), respectively. We also qualitatively discuss the 
additional steps that are required to interpret and use them in a massive theory. 



3.1 Coupling 

The canonical definition of the coupling is 



(6) 



In the massive theory, we choose 



.9a 



64tt 5 p+<5 (5) (pi - p 2 



P3j 



exp ( 9— ) (^3|M a (0 A ,m,A)|0i)| 



P2=P3 



64tt 5 p+<5 (5) (Pi -P2-P3)} {<h<h\V(g A ,m,A)\<f>x)\ 



(7) 



which differs from the definition in the massless theory by the factor exp 

(9t£)- This choice of 

coupling cancels a cutoff factor due to the presence of the mass and allows us to closely follow the 
formalism developed in the massless theory. In particular, the expressions for the matrix elements 
of VQ£,(m, A) and V^(m) presented below have the same form as those derived in Ref. Q. 
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3.2 Cutoff-Dependent Contributions to V^ r \m, A) 

Momentum conservation implies that any matrix element of V^ r \m, A) contains a sum of terms, 
each with a unique product of momentum-conserving delta functions. Assuming that approximate 
transverse locality is maintained, the coefficient of each product of delta functions can be written 
as an expansion in powers of transverse momenta. In massive (j) 3 theory, we can also make a 
generalized expansion in powers and logarithms of m. The scale dependence of any term in this 
expansion has the form 



A' 



-2jV lD 



m 
log x 



P± 
A 



(8) 



where A^; nt is the total number of particles in the final and initial states that participate in the 
interaction. Also a, (3, and 7 are non-negative integers. For simplicity we display one component 
of transverse momentum, however, the general form includes a product of all transverse com- 
ponents from all particles. In principle, the introduction of a particle mass allows any function of 
^ to appear. However, to O(g^) the only extra scale dependence comes in the form (^)° [log 7^] • 
If (3 = and 

6 - 2iV int _ a - 7 = 0, (9) 

the term is independent of the cutoff and is referred to as a "cutoff-independent" contribution. 
These contributions are discussed in the next subsection. 

The expression for a matrix element of Vq^ju, A) is derived from Eq. (|j): 



(F\V^(m,A)\I) 



(F\5V^ r \m, A, A')|/> - £ B r ^ s (F\V^ s \rn, A)\I 



(10) 



A terms 



"A terms" means the terms in the momentum and mass expansion that contain A' are to be removed 
from the expression in brackets. In terms that depend on positive powers of A', we do this by letting 
A' — > 0, and in terms that depend on negative powers of A', we let A' — > 00. 



3.3 Cutoff-Independent Contributions to V^ r \m, A) 

Considering the condition in Eq. @, only two-point and three-point interactions can have cutoff- 

(3) 

independent contributions. The lowest-order cutoff- independent three-point interaction is Vqj (m) 

(2) 

and has not been explicitly computed in the massless or massive theories. However, VX ^rri) is the 
lowest-order cutoff-independent two-point interaction and must be calculated before anything is 
calculated to third order. 

The matrix elements of (m) are divided into 2-point and 3-point contributions, and are 
given by the expression 



(F\V&\m)\I) 



1 



B r ,2 

1 



r+l 



(F\SV^ r+2 \m, A, A') I J) - J2 B r+2-s,s( F \V { 



r+2-s) 



m)\I) 



s=3 



+ 



B 



r.2 



(F\5V 



(r+2) 



r+l 



m,A,A')\I)-J2B r +2- s , s ( F \V {r+2 - s Xm)\I) 



3— point 

mPp9 term 
2— point 

m 2 term 



(11) 



Here, u m°p® term" and "m 2 term" means expand the term in brackets in powers of external trans- 
verse momenta and in powers and logs of m, and keep only the term that is proportional to m°p® 



4 



or m 2 , respectively.^ The removal of A and A' dependence is guaranteed by construction. 

Initially Eq. (|TT|) looks useless because Vqj (m) depends on V^j (m) [which is inside an 
integral in dV^ r+2 ^], suggesting the theory must be solved to all orders simultaneously. However, 
contributions to the reduced interaction from three-point interactions can only appear at odd orders, 
and contributions from two-point interactions can appear only at even orders. Thus, in the massless 
theory, this apparent problem does not manifest itself because there are no cutoff-independent 
two-point interactions. In the massive theory, although there are cutoff-independent two-point 

(2) (3) 

interactions, it is possible to solve for Vqj (to) and Vq T (m) simultaneously, without considering 
higher orders. This even-order /odd-order solution pattern can be extended to all orders. 

Including self-energy contributions, the theory we want to describe contains particles of mass 
to. We can simplify the problem by using this fact instead of using Eq. ( |TT1) to solve for the even- 
order VQj\m)'s. We do this by forcing the completely disconnected parts of the forward T-matrix 
elements to be zero. (This part of a T-matrix element contains initial and final states that have 
the same number of particles, n, and n momentum-conserving delta functions.) This fixes the 
even-order V^j\m)'s since they only involve interactions on single particle lines. This allows us to 

(2) (3) 

calculate V^j (to) independently of V^j (to). This extra condition can be used to fix all even-order 

4 Results 

The coupling in this theory runs at third order. We can compare the coupling at two different 
scales, A and A': 

oo 

We can determine how the coupling runs at third order by solving for Cs(to,A, A') (which is pro- 
portional to the matrix element ((/>2(/>3|<5W 3 \m, A, A') |p 2=P3 ). Figure [j] shows how C3(to,A, A') 
depends on the mass. The running of the coupling is exponentially damped as the mass grows since 
the cutoff inhibits production of intermediate particles. The difference between the values of the 
running coupling at two different scales increases as the two scales are separated. This is shown by 
the larger magnitude of C^{m, A, A') as the separation between A and A' grows. 

(3) 

Determining V^j (to) requires a fifth-order calculation and is not attempted. However, calculat- 

(3) 

ing the matrix element ((^^sl^cj^m, A) gives the relative sizes of the non-canonical interactions 
and the canonical interaction. Their relative magnitudes are similar to those in Ref. ffl], suggesting 
that an expansion of the reduced interaction in powers of the running coupling is valid through 
third order. 

Figure |2| shows how the non-canonical part of the matrix element of the invariant-mass operator 
for the interaction 0i — > 4>2<fi3 depends on the magnitude of the transverse momentum in the 
center-of-momentum frame. Increasing the transverse momentum in the center-of-momentum frame 
increases the free mass of the system. 

This work was partially supported by National Science Foundation grant PHY-9800964. 



The 2-point contribution is independent of p± because of cluster decomposition, transverse rotational 
invariance, and boost invariance. 
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Figure 1: The third-order coefficient of the 
running coupling function of the parti- 
cle mass. Curves for various upper cutoffs 
with fixed lower cutoff show the coupling is 
exponentially damped with increasing mass. 



Figure 2: The matrix element of the non- 
canonical part of the invariant-mass opera- 
tor for (pi — > 0203 versus the magnitude of 
the transverse momentum in the center-of- 
momentum frame, y is the longitudinal mo- 
mentum fraction carried by particle 2. 
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